Ultrasensitive iV-photon interferometric autocorrelator 
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We demonstrate a novel method to measure the iVth-order (N—l, 2, 3, 4) interferometric auto- 
correlation with high sensitivity and temporal resolution. It is based on the combination of linear 
absorption and nonlinear detection in a superconducting nanodetector, providing much higher ef- 
ficiency than methods based on all-optical nonlinearities. Its temporal resolution is only limited 
by the quasi-particle energy relaxation time, which is directly measured to be in the 20 ps range 
for the NbN films used in this work. We present a general model of interferometric autocorrela- 
tion with these nonlinear detectors and discuss the comparison with other approaches and possible 
improvements. 



The temporal correlation functions of various orders 
are of fundamental importance in the classical and quan- 
tum description of optical fields. The first-order (field) 
autocorrelation function describes temporal coherence 
and therefore spectral linewidth, the second-order (in- 
tensity) autocorrelation is used to measure the tem- 
poral properties of pulsed sources and to distinguish 
quantum and classical fields, while the measurement of 
higher-order autocorrelation is more sensitive to coher- 
ence features (e.g. photon bunching) of the light field [1] 
and can be used to determine the asymmetry of light 
pulses [2J. While the first-order autocorrelation func- 
tion is easily measured using an interferometer and a 
linear detector [3] , the measurement of higher-order cor- 
relation functions requires a process which is nonlin- 
ear in the intensity I(t). In interferometric autocor- 
relators, the normalized second-order correlation func- 
tion 3 (2) (t) = (I(t)I(t + T))/(I(t)} 2 is usually measured 
by using either second-harmonic generation (SHG) in 
a nonlinear crystal, followed by a linear detector 0], 
or two-photon absorption (TPA) in the detector itself 
[HE]. In both cases, the detector measures the square 
of the total intensity at the output of the interferome- 
ter (I t 2 ot (t)) oc (I 2 (tj) + 2(I(t)I(t + T d )) where r d is the 
delay between the two arms of the interferometer, to- 
gether with interference terms which are sensitive to the 
phase properties of the beam. While these approaches 
offer very high temporal resolution, since the related pro- 
cesses are nearly instantaneous, their sensitivity is lim- 
ited by the low nonlinear susceptibilities involved in the 
SHG or TPA process. Due to the even lower relevance of 
higher-order optical nonlinearities, they cannot be effec- 
tively applied to the measurement of the autocorrelations 
of order N>2. An alternative approach consists of com- 
bining linear optical detection with nonlinear processing 
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in the electrical read-out, e.g. in a correlation card, as in 
the Hanbury Brown and Twiss interferometer [7] . In this 
case, sensitive single-photon detectors can be used, how- 
ever the temporal resolution is limited > 100 ps by the 
jitter of the detector output and of the amplification and 
correlation electronics [8]. In addition, differently from 
interferometric autocorrelation, these approaches do not 
provide any information on the phase properties. 

In this work we report a novel approach to the mea- 
surement of the interferometric autocorrelation of order 
N > 2, which is based on the combination of linear 
absorption and nonlinear detection in a single device. 
The general principle consists of absorbing incident pho- 
tons in a linear absorber (i.e. a material where the ab- 
sorption probability per unit time is proportional to light 
intensity), which produces an output pulse only if two or 
more photons are absorbed within a certain time inter- 
val in the fs or ps range. We show the implementation of 
this concept in a superconducting nanodetector, where 
the nonlinearity is widely tunable by varying the bias 
current. We directly measure the temporal dynamics of 
the nonlinearity in the ps range and we attribute it to 
the relaxation dynamics of photo-created quasi-particles 
(QP). We show its application to the measurement of up 
to fourth-order interferometric autocorrelation, observ- 
ing an extremely high sensitivity related to the linear 
absorption process and to the low detector noise. 

The superconducting nanodetector consists of a 
nanoscale constriction in a superconducting wire (see 
FigjlJ. The device used in this work is based on a 
4.3 nm-thick NbN film (critical temperature T c =10.2 K) 
sputtered on GaAs substrate and has a constriction size 
of about 150 nm, patterned by electron beam lithogra- 
phy and reactive-ion etching. The nanodetector is bi- 
ased with a current Ih smaller than the superconducting 
critical current I c . Similarly to nanowire superconduct- 
ing single-photon detectors [TU] , the absorption of one or 
more photons produces a non-equilibrium population of 
QP in the nanodetector's active region, locally suppress- 
ing the superconductivity and increasing the probability 
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of vortex crossing, which can result in the transition to 
the normal state [11] [12] . The detection probability is a 
strong function of the absorbed energy (i.e. to the num- 
ber of photons) and of the bias current Jb, so that the 
nanodetector can be set to respond to > N photons by 
choosing lb 02]. A full tomographic characterization of 
the multi-photon response of the nanodetector under il- 
lumination with short pulses has been reported in Ref. 




FIG. 1: Scanning electron microscope image of a NbN super- 
conducting nanodetector. 

In this study we focus on the temporal characteristics 
of the nanodetector's multi-photon response and on its 
application as an interferometric autocorrelator. As we 
focus here on the application to the characterization of 
classical light sources, we describe the detection process 
as a sequence of single-photon absorption events using 
the semi-classical photo-detection theory. Assuming an 
incident light pulse with cycle- average intensity Ii n (t) on 
the nanodetector, the probability of creating a hot-spot 
within a time interval [t + dt] is equal to £ • S ■ Ii n (t) ■ dt, 
where C is related to the absorptance 7y a b s by £ = ij^bs/hv 
and S is the active area. In the two-photon regime, for 
example, the click probability P c i; c k (assumed -C 1) is 
equal to the probability that two photons are absorbed, 
weighted by a function 77(712) of the time difference t\% 
between absorption of the first and the second photon at 
times t\ and t 2 respectively, 

+00 +00 

Pciick = J J V(h-tl)CSIin{tl)CSIu,(t2)-dtfdt 2 

— OO — OO 

+ OO +OO 

= ( 2 S 2 J 17(7-12) J 4x(i)/in(i+7l2)-dt-d7l2, (1) 

— oo — oo 

The t7(ti2) , named nonlinear response function (NRF) 
in the following, depends on the QP dynamics in the su- 
perconductor and determines the intrinsic response time 
tnd of the nanodetector. The 77(712) is expected to decay 
from a maximum value of 77(0) (which depends on 7b) to 
for T12 3> tnd- In general, the detection probability in 



the iV-photon regime (N > 2) is, 

+00 +00 

Pciick = ( N S N J ri N (Ti N ,T 2 N, ■ ■ ■ ,Tn-1,n) J Iin(t) 
— 00 —00 
hn{t + njvKin(* + T 2N) ' ' ' ^in(* + T N -1,n) 

■dt ■ dr 1N ■ dr 2N ■ ■ • drjv-i.jV (2) 

Where t,-jv — tjf — U denotes the difference between 
the absorption times of the ith- and the last photon. 

When the incident pulse width is much larger than 
r ND , Eq.H is approximated as P c iick oc I^(t) ■ dt. 
When a nanodetector is placed at the output of a Michel- 
son interferometer, the input intensity reads I- m (t) oc 
\E(t) + E(t + T d )| 2 and the P c iick in iV-photon regime is 
proportional to the iVth-order interferometric autocorre- 
lations given by P c iick(Td) oc f+™ [\E(t) + E(t + r d )| 2 ] • 
dt. In the two-photon regime, by filtering out the inter- 
ference terms we obtain P c n c k( T d) (I 2 {t)) + 2(I(t)I(t + 
7d)), and therefore the g^(Td), similarly to the SHG- 
and TPA-based autocorrelators. 

The NRF can be measured by probing the autocorre- 
lator with short pulses. For pulse duration much shorter 
than tnd, the response of the nanodetector placed at 
the output of a Michelson interferometer can be found 
from Eq.Q. In the two-photon regime for example, af- 
ter filtering out the interference terms, Eq.Q becomes 
Pciick(Td)cx 77(0) + 77(r d ) + /(T d ) H3|. The function /(r d ) 
depends on the degree of first-order coherence of the in- 
put light and is different from zero only for delays shorter 
than the coherence time. For longer delays, the P c iick(Td), 
normalized by its value at 3> tndi is expected to vary 
as P c iick(7d)/Pciick(oo) = l + r;(r d )/77(0), so that the NRF 
and therefore the autocorrelator's timing resolution can 
be measured. 

The NRF was first measured by sending 1.6 ps pulses 
from an optical parametric oscillator (OPO) at A — 
1.13/im into a fiber-based Michelson interferometer and 
then to a nanodetector held at a temperature of 1.2 K 
using a lensed fiber producing a spot with an e~ 2 di- 
ameter of 5 /im. The delay r d in the interferometer is 
controlled by a motorized delay line (coarse control) and 
a fiber stretcher (fine control). The I c of the device was 
about 26 /iA. During the measurement, the nanodetec- 
tor was set in different photon regimes by choosing a 
proper 1^. As shown in Fig. [2|a), the count rate (CR) 
was measured as a function of the light power at differ- 
ent Jb- The solid lines with slopes of 1.04, 2.06, 3.06 
and 3.99 are the fits to the measured data in log-log 
scale, in the power ranges where the one-, two-, three- 
and four-photon response was dominant [13]. We chose 
one point at each of the four photon regimes, and mea- 
sured the CR as a function of r d , normalized by its values 
at long delays, as shown in Fig. [5Jb). The data points 
near r d = 0, where the measured autocorrelation is sen- 
sitive to the first-order coherence, are not shown in the 
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plot and were not considered in the fit since they intro- 
duce additional fitting error [15] ■ At p, = 18.0 fiA, the 
CR is independent of the delay, apart from the short- 
delay interference fringes, as expected, since the detector 
is working in the linear regime. When the p was low- 
ered to 12.5 ^A, 9.7 [iA and 8.4 ^A (corresponding to 
the two-, three- and four-photon regime, respectively), a 
maximum is observed at zero delay. Since the width of 
these peaks is much wider than the OPO pulse width, the 
measurement probes the intrinsic response of the detec- 
tor. In particular, the two-photon normalized response 
in Fig. [2jb) directly provides the 1 + T](Td)/f](0) depen- 
dence. The rj{r) is determined by the thermalization and 
relaxation processes of the photo-created QP and by the 
functional dependence of the P c ii c k on the QP concentra- 
tion jTTJ [T2] . Indeed, the QP population first grows as 
the electron population thermalizes via electron-electron 
scattering, in a timescale of few ps, then decays due to 
electron-photon interaction and phonon escaping to the 
substrate [TBJ. This QP decay is expected to determine 
the tnd in the TV-photon regime, since the p is cho- 
sen so that Pciick is high only when all QP produced by 
the iV photons are present at the same time. As P c iick 
is a strong function of the QP number, involving many 
microscopic parameters, a fit of the measured tj(t) us- 
ing a microscopic model would not be reliable. Instead, 
we introduce an empirical Gaussian NRF defined as 
7 7( T i2) = ??(0)exp[— (ti2/tnd) 2 ]- The fits to ten measured 
two-photon autocorrelation traces provide a tnd value of 
20.4 ± 0.8 ps. In order to fit the N > 2 traces we fur- 
ther assume that the multi-photon response factorizes as 

Vn(tin,t2n, ■ ■ ■ ,tn-i,n) = v{ t in)v( t 2n) ■ ■ -v( t n-i,n)- 
This is reasonable if one assumes that rj jy has an approx- 
imately exponential dependence on the total QP concen- 
tration after absorption of the TVth-photon, as suggested 
by the vortex-assisted photo-detection model jTTJ IT2"] . 
and that the QP relaxation time does not depend on QP 
concentration. Using the tnd value extracted from the 
two-photon autocorrelation traces as described above, 
Pciick (Td) was calculated for the three- and four-photon 
regimes from Eq.([2| without additional fitting parame- 
ters and shows excellent agreement with the experiment 
(Fig. [2jb)), which provides a strong experimental sup- 
port to our model. The value of tnd measured here at 
1.2K is comparable but shorter than the hot-electron en- 
ergy relaxation time in NbN microbridge measured by 
electro-optic sampling [T7] and by Terahertz spectroscopy 
[IE] - this difference may be related to the difference in 
film thickness and film properties. 

With the knowledge of the temporal resolution, auto- 
correlation experiments have been performed on pulses 
generated by a gain-switched 1.3^m diode laser with 
10MHz repetition rate and about 70ps pulse width. The 
one-, two-, three- and four-photon regimes were first 
found by choosing different p. At each p, CR were 
recorded as a function of in the TV-photon (TV=1, 2, 
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FIG. 2: (a) CR as a function of the light power for differ- 
ent 7b. The fitting lines with slopes of 1.04, 2.06, 3.06 and 
3.99 indicate the one-, two-, three- and four-photon regimes 
respectively, (b) Normalized CR (circles) as a function of ra, 
measured when the detector was working at the four chosen 
points (marked by large circles in (a)). The solid lines are 
calculation results (see text). 



3, 4) regime as shown in Fig. J^a-d), respectively. The 
fringe contrast ratio is observed to increase with TV, in 
good agreement with the theoretical values of 2, 8, 32 and 
128 for TV=1, 2, 3 and 4, respectively. In the one-photon 
regime (Fig. [3]^a)), the normalized first-order autocor- 
relation function g^'(Td) was calculated from the visi- 
bility of the interference fringes [5]. For TV > 1, higher- 
order intensity autocorrelation traces were obtained by 
applying a low-pass filter to the interferograms in Fig. 
|3^b-d). In order to explain the experimental results, we 
assumed that the incident light was a Gaussian pulse 
with a linear chirp. The electric field is modeled as 
E(t) = P exp[-41n(2)(l + iA)t 2 /r£], where E is the 
field amplitude, A is the linear chirp parameter [19] and 
r p is the full width at half maximum (FWHM) of the 
pulse. By fitting the measured |g^( T d)| hi Fig. [3^a) 
together with the second-order intensity autocorrelation 
(low-pass trace) in Fig. [3^b) , A and t p were determined 
to be 5.3 and 70.6 ps. The third and fourth-order inter- 
fcrometric autocorrelations were then calculated based 
on Eq.Q without additional fitting parameters. Their 
low-pass traces show a good agreement with the experi- 
ments, considering the very simplified assumption for the 
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FIG. 3: (a-d) Normalized CR as a function of Td in the different regimes at 7b of 18.0 fiA, 14.0/iA,11.5/^A and 9.6^A respectively. 
The |<? (Td)| (red line in (a)) and low-pass curves (red lines in (b-d)) agree well with the calculations (blue lines). Upper-left 
insets: CR as a function of incident power in log-log scale. The fitting slopes of 1.01, 1.96, 3.08 and 4.04 indicate the one-, 
two-, three- and four-photon regimes, respectively. The green circles indicate the points chosen for the measurements in the 
main panel. Upper-right insets: Expanded view of the interference fringes at small delays, plotted as a function of the voltage 
Vfs applied to the fiber stretcher. Lower-right inset in (d): FWHM of the fringes normalized by their period extracted from 
the upper-right insets in (a-d) as a function of N, showing a l/\/iV dependence. 



chirp. An enlarged view of each interferogram is shown 
as upper-right inset in each panel of Fig. [3j showing a 
clear narrowing of the fringes for increasing TV. As shown 
in the lower- right inset of Fig. |3fd), the FWHM of the 
fringes, normalized by their period, scales as approxi- 
mately 1/VN, which is a characteristic of multi-photon 
interferometry |20] and further confirms our conclusions. 

A comparison of the sensitivity of our nanodetector to 
the detection schemes used in conventional autocorrela- 
tors provides useful insight. For input pulses longer than 
t nd , Eq.|2| can be written as -P c iick = CND P? n (t)dt, 
where p n = I- m S is the incident power and Cnd = 
%/ 7r7 7abs 7 7(0) T ND/(/iz y ) 2 represents a nonlinear response ef- 
ficiency. A similar expression for P c n c k can be derived for 
nonlinear detectors based on SHG or TPA, where the cor- 
responding Cshg and Ctpa values are related to the SHG 
nonlinear conversion efficiency and to the TPA coefficient 
/?2 , respectively. In all three cases, the sensitivity can be 
defined by imposing that the CR is equal to the dark 
count rate Pdark in the detector. For a periodic sequence 
of square pulses with repetition rate / rep the condition 

/rep-Pclick > -Rdark leads to PpkPw > Pdark/C, where P pk 

and P av are the peak and average power, respectively. 
The two-photon autocorrelation trace in Fig. [3^b) was 
taken at P p kPiv = 5.6 x 1CP 17 W 2 , about seven orders of 
magnitude lower than the minimum reported PpkPiv for 
TPA using a GaAs PM tube [ST] and about two orders 
of magnitude lower than the lowest PpkPw reported in 
the second-order interferometric correlations using SHG 
[22] . Using the measured values of rj a hs = 1-5 x 10~ 4 , 
r](0) « 0.5 [2] and Pdark = 1Hz, we derive a sensitiv- 
ity of Pp^ in P a ™ in « 5.8 x 1CT 20 W 2 , corresponding to -4 
photons per pulse in our experiment. In higher-photon 
regime the advantage of using linear absorption is even 
larger. Indeed, the three-photon autocorrelation shown 
in Fig. [3^c) was performed at an average power of about 
InW, corresponding to P 2 k P av « 2.0 x 10" 21 W 3 , an im- 
provement of about twenty-one orders of magnitude over 



Ref. [13] ■ A^-photon interferometric autocorrelation for 
N > 3 has not been reported before to the best of our 
knowledge. 

The very high nonlinear response of the nanodetector, 
can be directly traced to the finite size and time dura- 
tion of the hot-spot created by the real absorption of 
one photon, as compared to virtual transitions involved 
in TPA. Indeed, a two-photon detection is triggered if 
the second photon is absorbed within the volume and 
time duration of the hot-spot created by the first pho- 
ton. This shows that a compromise exists between Cnd 
and tnd: for the nanodetector the tnd is determined by 
the QP relaxation time, while in TPA it is related to the 
lifetime of the virtual states associated to the TPA tran- 
sition, of the order of fs [5J. A similar compromise exists 
in SHG-based autocorrelators, where higher conversion 
efficiency requires a longer SHG crystal translating into 
a smaller phase-matching bandwidth and lower tempo- 
ral resolution [24] [25]. Interestingly, in an absorption 
thickness of 4.3 nm the nanodetector reaches a nonlinear 
efficiency C«1.7x 10 19 W J , several orders of magni- 
tude higher than /im-thick TPA absorbers and compara- 
ble to cm-long SHG crystals, showing the giant effective 
nonlinearity achieved through the combination of linear 
absorption and nonlinear detection. We note that a key 
advantage as compared to SHG-based autocorrelators is 
the wide optical bandwidth, limited only by the require- 
ment to operate in the desired A^-photon regime, which 
can be easily adjusted by varying the p,- This in principle 
enables the measurement of interferometric autocorrela- 
tion from the visible to the mid-infrared. 

The sensitivity of our autocorrelator is presently lim- 
ited by the low ?7 a bs, related to the spatial mismatch be- 
tween the incoming beam and the nanodetector's active 
area and to the small thickness of the NbN film. By 
focusing the beam with a high-numerical aperture lens, 
a much higher absorptance 7y a bs ~ 10 -2 can be achieved 
[26) . The integration of a plasmonic antenna and of a bot- 
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torn reflector [37] could increase ?7 a b s to the 10 1 range, 
leading to P^ in P a ™ n « 10~ 25 W 2 for the two-photon auto- 
correlation. On the other hand, increasing the detector 
area (as done in meander nanowire detectors [2S]) re- 
sults in f^bs oc L and 77(0) oc l^s/L (L is the nanowire 
length, ^hs is the hot-spot length) so that Cnd scales lin- 
early with 7/abs- The nanoscale nature of the detector is 
therefore crucial to reaching the ultimate sensitivity. Fi- 
nally, we note that the 20ps temporal resolution in our 
experiments is limited to the QP relaxation time in NbN 
films and could be much improved using different super- 
conducting materials, such as high-T c Y-Ba-Cu-0 films, 
where relaxation times ~1 ps were observed [29], opening 
the way to the characterization of TVth-order correlation 
functions in few-ps range with unprecedented sensitivity. 
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Supplemental Material: which is set in the two-photon regime and placed at the 

output of the interferometer, is found using Eq. (1) and 
In the general case, the response of the nanodetector, reads, 



Pciick(r d ) cx 2 J V (t 12 )G (2 Ht 12 ) • dr 12 + J V (t 12 )G (2 Ht 12 + r d ) • dr 12 + J V (t 12 )G (2) (t 12 - r d ) • dr 12 

— OO —CO 

+ CO +CO 

J r/(r 12 ) J E*(t)E(t + r 12 )E(t + r d )E*(t + r 12 + r d ) • dt ■ dr 12 

— oo — OO 

+oo 

+ y i /?(ri 2 )Fi(ri2,cjr d ) • dri2 + J ^12)^2(^2, 2ur d ) ■ dr 12 



+2Rc 



+00 



(SI) 



The first term on the right side of Eq. (SI) is a con- 
stant with respect to r d . The second and the third 
terms represent the convolution of the NRF with the 
second-order correlation function defined as G^(t) = 
f*™I(t)I(t + t) ■ dt, while the fourth term is sensitive to 
the phase-coherence of the input beam and gives rise to 
an additional peak in the filtered autocorrelation signal 



r 



at small delays. The last two terms on the right side 
of Eq. (SI) are interference terms which vary with tor d 
or 2oor d . Since they average to zero they can be easily 
filtered out using a low-pass filter. 

In the case where the pulse is much shorter than tnd> 
Eq. (SI) is written as, 



click 



(r d ) CX 7?(0) + 7?(T d ) + - 



v(0)\G^(r d )\ 



interference terms 



(S2) 



J GW(t 12 ) ■ dr 12 



Where G^\r d ) is the first-order correlation function 
defined as G^(r d ) =f+™E*(t)E(t + r d ) • dt. The third 
term on the right side of Eq. (S2), which comes from 
the fourth term in Eq. (SI), depends on the first-order 
coherence properties of the input beam, but, differently 
from the interference terms, does not disappear upon fil- 
tering and produces a peak in the filtered autocorrelation 
at zero delay. 

In order to confirm this prediction, we studied the 
short delay window in the experiment where the OPO 
short pulses (1.6ps, 1.13/im) were used to characterize the 
NRF of the nanodetector in the two photon regime. To 
this aim, the delay r d was controlled by a motorized delay 
line (coarse control) and a fiber stretcher (FS) (fine con- 
trol). A driving voltage Vfs (zigzag wave, ±5V, 20mHz) 
was applied on the fiber stretcher to change the optical 
path difference between the two arms. At each coarse 
delay, the counter recorded the count rate as a function 
of the Vfs- Due to the short pulse width of OPO, the 
interference fringes were observed at small coarse delays 



and an additional peak appeared in a small window near 
zero delay (blue points as shown in Fig. 4(a), also shown 
in an enlarged view in the inset), on top of the broad 
peak related to the detector NRF, as observed in Fig. 
2(b) of the main text and replotted as green circles in 
Fig. 4(a). A theoretical comparison based on Eq. (SI) 
is shown in Fig. 4(b). It considers a 20.4ps intrinsic re- 
sponse time of the nanodetector, and a probe pulse of 
1.6ps duration and 1.13/xm wavelength as used in the 
experiments, showing a good agreement with the mea- 
surements. We note that this narrow peak appears also 
in the filtered autocorrelation, shown as red lines in the 
insets of Fig. 4(a) and (b). However, as indicated by 
Eq. (S2), its width relates to the first-order correlation 
function, and therefore to the coherence time, and not 
to the pulse width. Indeed, we checked numerically that 
in a chirped pulse the narrow peak width scales with the 
chirp parameter. This confirms that tnd represents the 
limit to the temporal resolution for the measurement of 
the second-order correlation function. 
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FIG. 4: (a) The measured count rate as a function of the time delay, normalized by its value at long delays, (b) Simulation 
based on Eq. (SI). An enlarged view at small delays is shown in the insets. 



